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Abstract: We present the analysis of the signature 2j + $t(+X) via squark-squark production 
and direct decay into the lightest neutralino, pp —> qq — >■ jjXiXi(+X), in next-to- leading order 
QCD within the framework of the minimal supersymmetric standard model. In our approximation 
the produced squarks are treated on shell. Thus, the calculation of production and decay factor- 
izes. In this way, we provide a consistent, fully differential calculation of NLO QCD factorizable 
corrections to the given processes. Clustering final states into partonic jets, we investigate the 
experimental inclusive signature 2j + fx for several benchmark scenarios. We compare resulting 
differential distributions with leading-order approximations rescaled by a flat K-factor and examine 
a possible impact for cut-and-count searches for supersymmetry at the LHC. 
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1. Introduction 

Supersymmetry (SUSY) [1] is one of the most appealing scenarios for physics beyond the Standard 
Model (SM), which the ongoing experiments at the Large Hadron Collider (LHC) are searching for. 
SUSY, predicting new states with masses at the TeV-scale or below, provides an elegant solution 
to the hierarchy problem, and gauge coupling unification can be achieved naturally. In addition, 
many realizations of SUSY, particularly the Minimal Supersymmetric Standard Model (MSSM), 
provide a viable dark matter candidate assuming the lightest SUSY particle (LSP) to be stable due 
to R-parity conservation. Furthermore, the MSSM is in accordance with the measured values of 
the muon anomalous magnetic moment and of electroweak precision observables, and also with a 
light Higgs boson as indicated by the ATLAS and CMS experiments [2]. 
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Within the MSSM and assuming conserved R-parity, SUSY particles (sparticles) are produced 
in pairs, and related searches have been performed at LEP, the Tevatron, and the LHC using vari- 
ous final-state signatures. Due to their color charge, squark and gluino production typically gives 
the largest contribution to an inclusive SUSY cross section at a hadron collider like the LHC. As- 
suming the lightest neutralino Xi to be the LSP, produced squarks and gluinos eventually decay 
into the neutralino which leaves the detectors unobserved. This results in the general experimental 
signature of jets+missing energy, which is one of the signatures that has been searched for by the 
experiments at the LHC. Within the constrained MSSM (CMSSM) resulting limits can be used to 
exclude squarks and gluinos with masses below 1 — 1.4 TeV [3,4]. Exact exclusion limits, however, 
depend on the detailed structure of the underlying model parameters and become much weaker e.g. 
in parameter regions with compressed spectra where final-state jets do not pass the cuts applied in 
the experimental searches [5-7]. 

Precise and reliable theoretical predictions for squark and gluino production are necessary for 
several reasons: - to set accurate exclusion limits, - to possibly refine experimental search strategies 
in problematic parameter regions, and - in case of discovery, to determine the parameters of the 
underlying model [8]. The last point becomes more important as generic mass bounds for squarks 
and gluinos are pushed to higher values and many of the proposed sophisticated methods (see e.g. [9] 
for a review) for parameter determination might not be feasible with low signal statistics due to a 
rather heavy spectrum. 

Until now precision studies of sparticle production at the LHC have focussed on inclusive cross 
sections, without taking into account phase-space cuts that have to be applied in any experimen- 
tal analysis. Although these inclusive quantities are of fundamental interest, both for exclusion 
limits and for parameter determination, they are not directly observable in high-energy collider 
experiments. Furthermore, precise knowledge of distributions of the decay products might help to 
determine the fundamental parameters of the model [9, 10] or even permits the measurements of 
the spin of the new particles [11, 12] and thereby helps to discriminate SUSY models from possible 
other extensions of the SM with similar signatures [13,14]. Thus, a fully differential prediction 
including higher-order in all relevant stages of the process (eventually matched to a NLO parton 
shower) is desirable. 

In this paper we systematically study squark-squark production and the subsequent decay into 
the lightest neutralino at next-to-leading order (NLO) in QCD. Final state partons are clustered 
into jets and thus, we provide, for the first time at NLO, a fully differential description of the 
physical signature 2j + ]£t(+X) via on-shell squark-squark production and decay. In principle, 
our calculation does not depend on the hierarchy between the squarks and the gluino. However, in 
our numerical evaluation we only consider benchmark points where the mass trig of all light flavour 
squarks is smaller then the gluino mass m g > niq; otherwise the decay of a squark into a gluino 
and a quark would be dominant. Investigating the squark-squark channel should be understood 
as a first step towards a fully differential prediction for all sparticle production channels at NLO. 
It is, however, also of practical importance, since from recent searches at the LHC mass bounds 
for squarks and gluinos are generically pushed to higher values and here squark-squark production 
(initiated from valence- quarks) yields the dominant channel [15]. 

First leading order (LO) cross section predictions for squark and gluino production processes 
were already made many years ago [16-20] and are reviewed in [21]. Also the calculation of NLO 
corrections in perturbative QCD has been performed quite some time ago [22-25]. These corrections 
can be large (10% — 70%, dependent on the process and the parameters) and have to be included 
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in any viable phenomenological study due to the otherwise enormous scale uncertainties (including 
NLO corrections the scale uncertainty on inclusive cross sections is typically reduced to an order 
of 20% — 30%). Besides the scale uncertainty, PDF uncertainties dominate the error of theoretical 
predictions of sparticle production processes. In a very recent publication [26] general guidelines 
for the systematic treatment of these errors are presented. 

Already in [24] differential distributions at NLO are presented for the produced squarks and 
gluinos. For these results, NLO correction factors (K-factors) arc usually assumed to be flat in 
phase space and used as a global multiplicative factor to the LO cross section. Both, LO cross 
sections and NLO K-factors can be calculated efficiently with the computer code Prospino 2 [27]. 
However, a systematic study of the differential behaviour of these K-factors has never been per- 
formed. Furthermore, in [24] (and Prospino 2) NLO corrections for squark-squark production are 
always summed over the various flavour and chirality combinations of the produced (light-flavour) 
squarks, whereas realistic physical observables do depend on the chiralities through the decay modes 
which are in general quite different. Hence, in this work we treat the individual squark chirality 
and flavour configurations independently. 

More recently also results beyond NLO in QCD were calculated, based on resummation tech- 
niques [28-35]. These corrections increase the inclusive cross section by about 2% — 10% and further 
reduce the scale uncertainty. Moreover, electroweak contributions can also give sizeable corrections. 
At leading order they were first calculated in [36,37] and at NLO in [38-44]. In detail, those cor- 
rections depend strongly on the model parameters and on the flavour/chiralities of the squarks. 

A similar amount of work has been put into the calculation of higher order corrections to 
decays of (coloured) sparticles, with focus mainly on the integrated decay widths and branching 
ratios. NLO QCD corrections to the decay of light squarks into neutralinos and charginos were 
first calculated in [45,46] and to heavy squarks also in [46] and in [47]. Corrections to the total 
decay width of light squarks are in general moderate (below 10 %) and can change sign dependent 
on the involved mass ratios. However, for very small mass splittings between the decaying squark 
and the neutralinos/charginos these corrections increase significantly. Higher-order corrections to 
the decay of top-squarks are in general sizeable, but depend strongly on the mixing in the heavy 
squark sector. Related to this mixing also decays into weak gauge bosons or Higgs bosons can 
become relevant [48,49], receiving large higher order corrections [50-52]. Decays of a gluino into 
a light squark and a quark at NLO QCD together with the decay of a light squark into a gluino 
and a quark have been calculated in [53]. Corresponding decays involving top squarks stops were 
presented in [54]. All these decays including their NLO QCD corrections have been implemented 
in the public computer programs SDECAY [55] and SUSY-HIT [56]. 

Besides NLO QCD, also NLO electroweak corrections to squark decays into neutralinos and 
charginos have been investigated in the literature [57,58] and can give sizeable contributions. These 
corrections often compensate those from QCD on the level of integrated decay widths, however, 
they depend strongly on the model parameters. Corresponding NLO electroweak corrections for 
third generation squark decays have been studied in [59-62]. 

As already mentioned before, most of the discussed studies of higher-order corrections focussed 
on inclusive observables or considered differential distributions in unphysical final states, like unsta- 
ble sparticles. Few studies were performed investigating invariant mass distributions of SM particles 
emitted from cascade chains including various higher-order corrections [63,64]. Finally, in [65,66] 
production of sparticles was studied at tree level matched to a parton shower including additional 
hard jets. In these works large deviations from the LO prediction with or without showering were 
found particularly in the high-p T tail for scenarios with compressed spectra. In this paper we 
go beyond these estimates and provide a fully differential description of production and decay of 
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squark-squark pairs at NLO. 



In the big picture of the complete calculation of NLO QCD corrections to pp — > 2j + 
also qq'*, gq and gg intermediate states can contribute to this signature. Already without system- 
atically including decays of the squarks, the calculation of NLO corrections to on-shell production 
of such pairs of coloured sparticles carries problems of double counting. Parts of NLO corrections 
to one final state can be identified as LO of another final state where the decay is already included. 
The standard solution, used to avoid this double counting problem, can not be straightforwardly 
extended to the calculation where decays (and off-shell effects) are included from the beginning. 
Moreover, complete NLO corrections to pp — > qq'xiXi do not only include factorizable contribu- 
tions, i.e., contributions that can be classified as corrections to the production or to the decays, but 
also non-factorizable contributions, where such a separation is not possible. In this paper we ana- 
lyze the factorizable NLO corrections to squark-squark production and decay, which are expected 
to yield the dominant part of the NLO contributions. Non-factorizable effects and off-shell contri- 
butions will be analyzed in a forthcoming publication, providing a consistent conceptual approach 
and evaluating their numerical effects. 

The outline of this paper is as follows. In section || the method of combining consistently 
production and decay at NLO in the narrow-width-approximation is described. In the following 
two sections the calculation of all required ingredients of this combination is explained, with respect 
to the squark production processes in section ^, and to the squark decays in || In section ^| we 
present our numerical result for representative benchmark points, and conclude with a summary in 
section O. 



2. Method 



We investigate the production of squark-squark pairs induced by proton-proton collisions, with 
subsequent decays of the squarks into the lightest neutralinos. Since we are interested in the exper- 
imental signature 2j + J£t(+X), all contributions from light- flavour squarks have to be included, 
hence the cross section is given by 

do = [ d<T (PP ~> Qi,aQj,b -> <HXi1jXi(+ x )) + do{pp -> q^ ta Q^,b ^ QiXiQjXii+X))] . 

a,b£{L,R} 

(2.1) 

Indices i,j denote the flavours of the (s)quarks and a,b their chiralities. At LO, the only partonic 
subprocesses that contribute to a given intermediate configuration qi^ a <jj,b or q* a q* b arise from quark 
and anti-quark pairs, respectively, qiqj -> -> qiXiljXi an d <Mj ~> Qi, a 9j,b ~> QiXiQjXi- 

For simplifying the notation, we will write qq' — > qq' — > qx\q'x\ whenever the specification 
of flavour and chiralities is not required . Moreover, we will perform the discussion without the 
charge-conjugate subprocesses; they are, however, included in the final results. 

In the considered class of processes, squarks appear as intermediate particles with mass nig and 
total decay width F^. In the limit Tq/rrig — > , the narrow width approximation (NWA), their 
resonating contributions in the squared amplitude can be approximated by the replacement 

1 7T 

9^5 9-^9 -> W 6 {P 2 ~ m d) ! ( 2 - 2 ) 

( P 2 - TO?) 2 + m? r? Tflq Yjj ^ «" V ' 

1 In this notation q = q' implies q = q', but not vice versa 
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for each squark with momentum p. 

At LO in NWA for scalar particles, the phase-space integration of the squared amplitude for 
the total cross section of the 2 — > 4 processes factorizes into a production and a decay part. At the 
partonic level, the LO cross section gets the following form, 

with the LO partonic production cross section v^qi^qqi and the LO branching ratios BR(°) for 



the squark decays into the lightest neutralino. A direct generalization of eq. (2.3) yields the cross 



section in a completely differential form, which can be written at the hadronic level as follows, 

^nwa^ -+ & vxWxl) = da^ r i drfl ? . ; -L drg^ . (2.4) 

<z 9' 

Therein, r~ ^ and T„ denote the LO total widths of the two squarks; dT.\ - and dTv 1 , - are 
the respective differential decay distributions boosted to the moving frames of g and <f . 

The other basic ingredient, da^J^--, , is the hadronic differential production cross section, ex- 
pressed in terms of the partonic cross section da^ q ,^~-, as a convolution 

d*£> - I drC qq .{r)da^^ ss ,{r) (2.5) 



99 V ; "'"gg' 

TO 



with the parton luminosity 



1 /"^ c?x r st\ 
Cqq'(r) = —— / — f q (x,fi F )f q > (-,10? + ( 9 « g ') , (2.6) 



where fi(x,/j,F) is the parton distribution function (PDF) at the scale /if of the quark i with 
momentum fraction x inside the proton, r denotes the ratio between the squared center-of-mass 
energies of the partonic and hadronic processes, r = s/S, and the kinematical production threshold 
corresponds to To = (m q + m^) 2 / S. 

The NWA cannot be extended to the complete set of NLO QCD corrections to pp —> qq 1 —> 
qXil'Xv Interactions between initial and final state quarks, for example, do not allow to split 
the process into on-shell production of squarks and subsequent decays. The subset of factorizable 
corrections, however, can be obtained in NWA, occurring as corrections to the production or to the 
decay processes. In the present article we focus on this class of corrections, which are expected to 
be the dominant ones. In an upcoming article we will investigate the non-factorizable corrections 
and their numerical influence. 

The differential cross section for on-shell production of squark-squark pairs and subsequent 
decays including the NLO factorizable corrections can be written as a formal expansion in a s , 



r (i) r (i) . 



<$ a (pp - m' - ti&m+x)) = ^ [d*^ w dvf^ (1 - -4-4) 



9 9' 



pp->gg g— >gx" g — >1 Xi pp->-qq g-^gx" 9 — >q x" 



pp^qq'(X) g->gx? 9'^9'X? 



(2.7) 



with the NLO contributions to cross section and widths da^,T^\ ... in obvious notation. The 



LO term in the first line of eq. (2_/7) gets a global correction factor from the NLO contribution to 
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the total widths; the second and third line involve the NLO corrections to the decay distributions 
and the production cross section 2 . 



In order to evaluate the terms contained in eq. (2.7), we produce, for all different combina- 
tions of light flavours and chiralities, weighted events for squark-squark production and squark 
decays. Production events for pp — > qq'(X) are generated in the laboratory frame. Decay events 
for q — > qxl(g) and q' — » q'xiid) are generated in the respective squark rest frame. Finally, 
pp —5- qq' XiXi(+^0 events are obtained by boosting the decay events from the squark rest frames, 
defined by the production events, into the laboratory frame. The weights of the pp — > qq' XiXii+X) 
events are obtained combining the different LO and NLO weights of production and decay according 
to eq. ( j2.7| ). Phenomenological results derived by these combinations are presented in section |^. 
The treatment of the various entries in eq. (2/7) is described in the following sections || and |j 



3. Squark-squark production 
3.1 LO squark squark production 

Amplitudes and cross sections for squark production depend on the flavours (indices and on 
the chiralities (indices a, b) of the squarks. Since we consider light-flavour squarks only, all quarks 
can be treated as massless. 



fJ 



1j,b 



Figure 1: Tree level Feynman diagrams for squark-squark production. 



If the two produced squarks are of the same flavour, the contributing Feynman diagrams corre- 
spond to t- and u-channel gluino exchange (figure [j]) . For squarks of the same chirality, the partonic 
cross section reads as follows, 



do 



(o) 

dt 



1 



f 



2/3 



9s 2 



(t — to?) 2 (u — m|) 2 (u — m?)(i — to?) 



(3.1) 



where s,t and u are the usual Mandelstam variables for 2 — > 2 processes. For different chiralities, 
m q = implies vanishing interference between the t- and u-channel diagrams, yielding 



da x ?l. lS . z. 2?ra 2 f st - (t - m 2 ~ za ){t - m? J -su - (u - m| <a )(u - m| 4 J 



dt 



9.s 2 



(£-m?) 2 



(u — TO?) 2 



(3.2) 



a' v g' 

If the two squarks are of different flavours, there is no it-channel exchange diagram; the partonic 
cross section for equal chiralities is hence given by 

- i - nr- s 

(3.3) 



do 



qiqj^qiaqja 



2nai 



dt 



9s 2 (t - ml) 2 ' 



and for different chiralities by 



do 



(o) 

dt 



27ra 2 (-t + mlJ(t-ml b )-st 



9s 2 



{t-ml) 2 



(3.4) 



2 An analogous treatment has been used, e.g., for the calculation of NLO corrections of top pair production and 
decay [67,68]. 
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Besides the dominating QCD contributions, there are also tree-level electroweak production chan- 
nels [36,43] with chargino and neutralino exchange, which can interfere with the QCD amplitude 
providing a contribution to the cross-section of 0(aa s ). In principle these terms can be numerically 
of similar importance as the NLO QCD 0{oP s ) corrections we are investigating. For the present 
study, the electroweak contributions are neglected. 

3.2 NLO squark squark production 

The NLO QCD corrections to squark-squark production have been known for many years [24] and 
an efficient public code (Prospino 2) is available for the calculation of total cross sections at NLO. 
However, in order to study systematically the 2j + ]St{+X) signature emerging from production 
of squark-squark pairs and subsequent decays into the lightest neutralino, also the complete dif- 
ferential cross section is necessary. To this purpose, we perform an independent (re)calculation 
of the NLO QCD corrections, where we treat the masses for q^, qn and all chirality and flavour 
configurations independently. In [24] different squark chiralities are treated as mass degenerate and 
NLO contribution are always summed over all chirality and flavour combinations. 

NLO calculations involve, in intermediate steps, infrared and collinear divergences. Since our 
calculation does not involve any diagrams with non-Abelian vertices, infrared singularities can be 
regularized by a gluon mass (A) . Collinear singularities, in analogy, can be regularized by a quark 
mass (m q ), that is kept at zero everywhere else in the calculation. The cancellation of these two 
kinds of singularities is obtained by summing the virtual loop contributions and the real gluon 
bremsstrahlung part, with subsequent mass factorization in combination with the choice of the 
parton densities. 

The complete NLO corrections to the differential cross section can be written symbolically in 
the following way, 

d<T ( V = dCT virtual+soft da ™u + da^ gg , a + da e ^^± . (3.5) 

pp-tqq'(+X) pp— >qq '(g) PP~ >qq Iff) PP~ ><J9 9 pp-tqq'qi'l <■ ' 

With <i<7 VIItua ^ s ? ft we denote the summed contributions from the renormalized virtual corrections 

pp^qq'{g) 

and soft gluon emission; dcr co11 \ corresponds to initial state collinear gluon radiation including 

o ' pp^rqq (g) * b ° 

the proper subtraction term for the collinear divergences; da^^--, g denotes the remaining hard 
gluon emission outside the soft and collinear phase space regions, da^J^'l^ is the contribution 
from real quark emission from additional quark-gluon initial states contributing at NLO. 

Technically, the calculation of the loop corrections and real gluon contributions is performed 
separately for every flavour and chirality combination, qiq^ — > qi a qjb, with the help of FeynArts [69] 
and FormCalc [70,71]. Appendix ^ shows a collection of the contributing Feynman diagrams. Loop 
integrals are numerically evaluated with LoopTools 2.7 [70]. 

3.2.1 Virtual corrections and real gluon radiation 

In the term da™ 1 ^^^ 1 the virtual and soft contributions are added at the parton level, according 
to 

dcr™^ ua l+ s ° ft = [ dr £w(r)da™ t " a i+ s ° ft (r), 

pp^qq'(g) I 11 V > qq'-yqq'(g) v n 

d °t'*tZ)^ = da^%,(m q A) + dalf^ s , (g) (m q ,X,AE). (3.6) 
The fictitious gluon mass A for infrared regularization cancels in the sum of da™^—, and da^^^-^i^ . 

At NLO, UV finiteness requires renormalization by inclusion of appropriate counterterms, which 
can be found explicitly in [39] . All mass and field renormalization constants are determined accord- 
ing to the on-shell scheme. The renormalization of the QCD coupling constant (dg s = g s SZ gs ) has 
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to be done in accordance with the scheme for a s in the PDFs, the MS scheme with five flavours; 
this corresponds to the renormalization constant [24] 



SZ a 



47T 



A 



00 



1 m 

q l0 S — 

3 fj,' 



log- 



1 

12 



9 



(3.7) 



with the UV divergence A = 2/e — 75 + log(47r) and the renormalization scale fi. /3q = 3 is the 
leading term of the f3 function for the QCD coupling in the MSSM. Applying the MS scheme we 
are forced to use dimensional regularization for the calculation. This breaks the supersymmetric 
Slavnov- Taylor identity that relates the qqg vertex function and the qqg vertex function at one-loop 
order. However, this identity can be restored (see [24,72]) by an extra finite shift of the g s coupling 
in the qqg vertex with respect to g s in the qqg vertex, 



fjs 



g s {l + SZg 



SZg s = SZg s 



3tt 



(3.8) 



The second term da 



oft 

qq' ->?<?' (9) 



in eq. (3.E) contains the contributions from real gluon emission 



integrated over the soft-gluon phase space with E g < AE. It is similar to the case of soft-photon 
emission [73,74], yielding a multiplicative correction factor to the LO cross section. In the case of 
gluons, however, the color structures are different for emission from t and u channel diagrams and 
hence the various bremsstrahlung integrals enter the cross section with different weights. Accord- 
ingly, we decompose the partonic LO cross section for qq' 
notation, 



qq' in the following way in obvious 



da 



(0) 



da 



(ut) 



^qq> 



n (ut) 

qq' 



qq' 



UO qq<^qq> 



(3.9) 



where the coefficients ci'~) for the individual channels can be easily read off from the LO cross 



sections in eqs. (3.1)— (3.4). Defining = 1 for incoming and e; = —1 for outgoing particles, 
the soft gluon contribution at partonic level can be written as follows, using the label assignment 
{<7,g',g,9'}o {1,2,3,4}, 



*<?«' ■ 



(3.10) 



The lij involve the bremsstrahlung integrals and the weight factors C~~, 



(tt,ut,uu) 



Explicit expressions 



are listed in eq. (B.l) of Appendix 



^pp-^lf'fg) depends on the quark mass m q owing to the initial-state collinear singularities. 
This dependence cancels by adding the real collinear radiation term da^^--, , n resulting from gluon 
emission into the hard collinear region and mass factorization for the PDFs via adding a proper 
subtraction term, 



J coll 

U(J pp^qq'(g) 



j coll-conc 
a "pp^qq'(g) 



da 



sub-pdf 

' pp->qq' ' 



(3.11) 



The collinear gluon emission into a narrow cone around the emitting particle yields the following 
contribution that corresponds to the results of [75] with the replacement aQ 2 q —> (4/3) a s , 



1 coll-conc 
a<J pp^qq'(g) 



1 dx 

X 



dz_ 

z 



' h I - I - ^(r, z) da^l? {g) (r, z) 



(3.12) 



The luminosity C c ° q , and the partonic cross section da coY[ " 
of Appendix [Br 



can be found in eqs. (BJ3) and (B.4) 
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The subtraction term for phase-space slicing, in accordance with the MS scheme, can be written 
in the following way, 



-2 f drC qql (r)da s ^ 9 ,(r), (3.13) 

where da^^^-i and c?<7™^~-, are defined in eqs. ( jgj) and of Appendix g. 

Finally, we have to add der har< J~~,, * from the residual hard gluon emission outside the collinear 

pp-tqq (.g) b 

region, which cancels the dependence on the slicing parameters S s for the soft and 5$ for the collinear 
region. 

3.2.2 Real quark radiation 

The 0(a s ) corrections to pp — > q~i ta q~j t b(-{-X) get also a contribution from the gluon-initiated sub- 
processes qtg — > qi, a Qj,bQj and qjg — > qi,a,qj,b<ji-, which also have to be included for a consistent 
treatment of NLO PDFs. Diagrams for these two subprocesses can be divided into resonant (fig- 
ure ||(a) and||(a)) and non-resonant (figure ||(b) and||(b)) diagrams, where in the resonant diagrams 
the intermediate gluino can be on-shell. This resonant production channel via on-shell gluinos cor- 
responds basically to LO production of a squark-gluino pair (with subsequent gluino decay). In 
contrast, the contribution with off-shell gluinos can be considered in different ways. In the more 
general context of production and decay of all coloured SUSY particles including off-shell effects 
everywhere, squark-squark and squark-gluino processes have to be treated together at NLO, and 
a separation can appear only as an intermediate organisatorical instrument. In the following we 
explain this strategy. 

In the case of different flavours i ^ j, there are two parton processes which, in our approach, 
provide NLO differential cross sections for real quark emission, given by 



d a qig^qi, a q jtb qj 



^qjg^qi,aqj,bqi 



\M nonres, % nonres. iM* TeSti ) , (3.14) 



l-Mnonrcsjl + 2Re(M nonre8 , i .M* esJ ) 



where overline represents the usual summing and averaging of external helicities and colours. M les ,i 
and nonres, i include respectively diagrams of figure ||(a) and figure ||(b), M ICS ,j and A^nonrcsj" 
include diagrams of figure |(a) and figure |(b). Here, we are excluding the squared contributions 

|A4res,i| and \M res, 7 

In an analogous way, for equal flavours i = j, we have 
da, 



'qig^qt 



\M nonres r + 2Re(M „o„resX* e J , (3.15 



where A^ rC s includes diagrams of figure ||(a) and figure ||(a), and A^ n0 nres includes diagrams of 
figure ||(b) and figure ||(b). Again, the term lA'lrcsl 2 is excluded. 

Such a splitting between resonant and and non-resonant contributions for off-shell gluinos has 
to be understood, as already mentioned above, as an organisatorical instrument. Only when the 
gluino is on-shell, resonant diagrams alone form a gauge-invariant subset. However, a consistent 
description is achieved when the prescription presented here is combined with relevant off-shell 
contributions. These will be studied in a forthcoming paper where we include the full contributions 
from the terms |7W ros | 2 , excluded here, together with off-shell effects of any intermediate coloured 
sparticle in all squark and gluino production channels. 
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(a) (b) 



Figure 2: Resonant (a) and non-resonant (b) diagrams contributing to qig — > qi, a qj,bqj- 




Figure 3: Resonant (a) and non-resonant (b) diagrams contributing to q^g —¥ qi, a qj,tqi- 

In the "Prospino scheme" [24,76] contributions from squared matrix elements with resonant 
diagrams are eliminated in a different way. A small non-physical width is used to regularize and 
subtract on-shell gluino-squark production with associated gluino decay. The "Prospino scheme" 
is well suited for calculations of production processes. It can not straightforwardly be extended to 
calculations where decays and/or off-shell effects are included in all channels. 



For the practical calculation of the real quark radiation contributions, one has to perform the 
phase space integration over the final state quark. The squared non-resonant terms in eq. ( 3.15 ) 
and eq. ( 3.14 ) lead to initial state collinear singularities. Again, these singular terms have to be 
subtracted since they are factorized and absorbed into the PDFs. Like in the case of gluon radiation, 
we divide the emission of a quark into a collinear and a non-collinear region (since no IR singularities 
occur, a separation into soft and hard quark emission is not required), 



i real-quark \ * 

a<J qia,qib<li/j ~ 2-^1 \ 



1 



j coil-quark . j noncoll-quark 
aa pp^.q ia q jb q k a(J pp^q ia q jb q k 



(3.16) 



The non-collinear contribution 



da 



noncoll-quark 

' pp^q ia q jb q k 



j noncoll-quark / \ ja 

dT Ajfc ( r ) dvuv^qiaqm* 



(t), 



(3.17) 



contains C 
r t 

C( 

J pp^tqia.qjbqk 



ijk C ° 11 quark (r) as given in eq. (B.9). The collinear emission together with the subtraction 



terms for the PDFs instead can be written as follows, 

1 rl dx f 1 dz 



da' 



coll-quark 



(5 ik + Sjk) 



, ^ . „™ -coll-quark/ x ,-coll-quark / x 

" T ' ' l ~ijk V ' ' / q%P — }Qin.Q4hQk \ ' ) ' 



tg^qiaqjbqk 



(3.18) 



with C ijk (T, x, z ) coll -q uark and da c ° l ^^ b - k (r, z) defined in eq. (^) and eq. (Q of Appendix | 
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4. Squark decay 



4.1 Squark decay at LO 



The LO decay width for a squark decaying into a neutralino and a quark, q — > qx®, depends on the 
flavour and chirality of the squark. For m q — the width can be written as follows, 

r(°) - " ( 



1 



x° 



f 2 

J a 



(4.1) 



The coupling constants f a can be expressed in terms of the isospin and the charge e q of the 
quark, together with the neutralino mixing matrix (Njk) including the electroweak mixing angle 
through % = sin(9w and cw = cos9w , 



f L =V2 e q N' + (ll L - e q s 2 w ) 



1 



CwSW 



N' 32 



f R = -V2\e q N' 



e q JV l'2 

c w J 



N' =c w N jl + s w N ]2 , NL = -swNj! + c w N. 



Vj2 ■ 



(4.2) 

(4.3) 
(4.4) 



For a scalar particle decaying in its rest frame there is no preferred direction, and hence the dif- 
ferential decay distribution is isotropic. For squark decays into neutralino and quark, the decay 
distribution is thus simply given by 



= — rl 0) . dcos6 
q^iXj 47]- q^qx" 



(4.5) 



with polar angle 9 and azimuth <j) referring to the quark momentum. 
4.2 NLO squark decay distribution 

The differential decay width for q — > qjft at NLO is obtained in analogy to the steps in section 



3.2 by adding the virtual loop corrections and the real gluon bremsstrahlung contribution from the 
soft, collinear, and hard non-collinear phase space regions, yielding the full NLO contribution in 
the form 



rfr 



jjHvirtual 
dL q^qx" 



j P soft 
ai -q^qX°(9) 



, jpcoll , jjnh 

+ dL q^qX°Ag) +dl q 



hard 



(4.6) 



The virtual corrections cTY u ' tu ?o for m. 



- 9 — correspond to the two vertex loop diagrams in figure|4](a) 

and the vertex counter term (indicated by the cross in figure |^(a)), which consists of the wave- 
function renormalization constants of the external quark and squark line. As for the production 
amplitudes, the renormalization constants are determined in the on-shell renormalization scheme. 
Details on the vertex counter term can be found in [39], and the explicit analytical expression is 
given m eq. ( PI of Appendix g 



qi.a 




qi, a 





q X ) 



q 




q 



fj 



9 

(a) (b) 
Figure 4: Loop and counterterm diagrams (a) and gluon radiation diagrams (b) for squark decays. 
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The term dT s ^^ can be calculated by the same strategy as for da g °, _^.gg, Q , yielding 



(0) 



q^qx 



(4.7) 



: ,i=i;i<i 



where the are defined in eq. flB.2| ). 

Collinear divergences now emerge from the final state. Making again use of the results of [75], 
the collinear emission of gluons with energy larger than AE into a cone with opening angle A9 
yields the contribution 



= drl°> _ 

q^-qXj 



(4.8) 



2a s 
3tt 



log 



2£ - 9 ,max 



•2 log (a.) 1-log 



op2 

g,max ' 



where <5 S = 2AE/niq, S$ = 1 — cos (A0) ~ A# 2 /2, and E qmax = — the maximum energy 

available for the quark in the squark rest frame. Gluons with 9 < A9 are recombined with the 
emitter quark into a quark with momentum p reC omb = Pq + Pg- In general, differential distributions 
in the quark momenta are dependent on the slicing parameter A9. However, this dependence will 
disappear once a jets algorithm that is much more inclusive in the recombination of quarks and 
gluons, is applied (see section ||). 

The contribution dr^ 1 ~ from real emission of hard gluons (E g > AE) at large angles (9 > AO) 

q qXj 

is evaluated by numerical integration of the squared matrix elements, obtained from the diagrams 
in figure 0(b). 



4.3 Total decay width 



The total squark decay width Iy' at LO is obtained by summing the partial decay widths of the 6 
different possible decay channels into neutralinos and charginos (assuming always nig > rriq). The 



partial decay widths into neutralinos are given directly by eq. (4.1). For charginos, the partial decay 
widths r 



(o) 



are also described by the formula (4.1), with the specification /r = and 



h = 



sw 



for 



u,c, 



h 



s w 



for 



d, 5, 



(4.9) 



for the coupling constants. Here, U, V are the mixing matrices in the chargino sector. 

For the total decay width at NLO, one has to calculate the NLO QCD corrections for each 
channel, which can be done analytically performing the full phase space integration over the three- 
particle final state with the radiated gluon and adding the loop contributions. The six partial decay 
widths contributing to r~ 0+1 ' at NLO can be expressed in terms of their respective LO result and a 
NLO form factor F® CD , which depends only on the mass ratios of the involved massive particles, 

4 a, 



-.(0+1) 
q^qx 

-.(0+1) 



,(0) 



-(0) 



1 



3 7T 



pQCD 









m g) . 



1 



i^lpQCD 
3 7T 









V 9 





(4.10) 



The result for F® CD in [46] is confirmed by our independent derivation. Details can be found in 
Appendix |c]. 



- 12 - 



5. Phenomenological evaluation 



Let us now turn to the numerical evaluation of of the process under consideration. In the following, 



we first specify in 5.1 the input parameters and benchmark scenarios we consider in our numerical 



evaluation. Then, in 5/2 a brief introduction to all considered distributions and observables is given. 
Finally, in 5.3, we present our numerical results in three different ways. Firstly, we compare our 
results for inclusive production of squark-squark pairs with Prospino 2 and investigate inclusive 
K-factors for different chirality and flavour combinations. Secondly, we present several differential 
distributions for various benchmark scenarios and center-of-mass energies. Thirdly, we investigate 
the impact of the considered higher-order corrections on total event rates and thus on cut-and-count 
searches for supersymmetry at the LHC. 

5.1 Input parameters 

Standard Model input parameters are chosen according to [21], 

M z =91.1876 GeV, M w = 80.399 GeV , G F = 1.16637 • 1CT 5 GeV , 
m^(m 6 ) = 4.2 GeV , m t = 173.3 GeV , m T = 1.777 GeV . (5.1) 

We use the PDF sets CTEQ6.6 [77] interfaced via the LHAPDF package [78] both for LO and 
NLO contributions. The strong coupling constant a^ s ((in) is also taken from this set of PDFs. 
Factorization scale (if and renormalization scale (ir are, if not stated otherwise, set to a common 
value, (i = (if = (J-R = Triq, with rriq being the average mass of all light-flavour squarks of a given 
benchmark point. 



For SUSY parameters, we refer to three different benchmark scenarios. First, we investigate 
the well studied CMSSM parameter point SPSla [79]. Although being practically excluded by 
recent searches at the LHC [3,4,80], this point still serves as a valuable benchmark to compare 
with numerous numerical results available in the literature. Second, we study the benchmark point 
CMSSM10.1.5 introduced in [81]. Due to its larger m\/2 parameter, compared to SPSla, squark 
and gluino masses are considerably larger, resulting in a generally reduced production cross section 
at the LHC. Not excluded yet, this parameter point can be tested in the near future. The overall 
spectrum, though shifted to larger masses, is very similar to the one of SPSla. Third, we consider 
a phenomenological benchmark point defined at the scale Q = 1 TeV. We follow the definitions 
of [81], where such a point sits on a line called pl9MSSMl. It can be parametrized by essentially 
one parameter, the gaugino mass parameter Mi. A unified parameter for the gluino and the light- 
generations sfermion soft masses M 3 = m n st /2nd gm is fixed to M 3 = m p st /2nd g on = 1.2 Mi on this 

line and we choose Mi = 300 GeV for our benchmark scenario pl9MSSMlA. All other masses and 
parameters as well as the soft masses for the first two generation right-handed sleptons are at a 
higher scale and irrelevant for our analysis. This benchmark point is chosen to study a particular 
parameter region with rather light squarks and gluinos which is not excluded by recent results from 
the LHC. Due to a small mass splitting between the Xi an d the light squarks (and gluino) resulting 
jets tend to be very soft and thus escape the experimental analyses. Particularly in such parameter 
regions precise theoretical prediction of the resulting SUSY signal including higher orders on the 
level of distributions seems to be necessary for a conclusive study. 

Parameters of the CMSSM benchmark scenarios are defined universally at the GUT scale and 
are shown in table |l| They act as boundary conditions for the renormalization group running of 
the soft-breaking parameters down to the SUSY scale MsuSY- This running is performed with 
the program SOFTSUSY [82] which also calculates physical on-shell parameters for all SUSY mass 
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benchmarkpoint 


m 


m 1/2 


A 


tan (3 


sign(/i) 


SPSla 


100 GeV 


250 GeV 


-100 GeV 


10 


+ 


10.1.5 


175 GeV 


700 GeV 





10 


+ 



Table 1: High energy input parameters for the two considered CMSSM scenarios. 



benchmarkpoint 


Mi 


M 2 


M 3 


Ai 


tan (3 


sign(/x) 


P19MSSM1A 


300 GeV 


2500 GeV 


360 GeV 





10 


+ 



Table 2: Low energy input parameters for the pl9MSSMlA scenario. The first two generation sfermion 

soft-masses (apart from the right handed sleptons) equal the gluino mass m -ist/2nd gen = M- A . All other 

'l/r 

parameters are at a higher scale m E , n = rng R ,p_ R = \i = Ma = M%. 

' h/R 



benchmarkpoint 


u L 


u R 


d L 


d.R 


g 


X? 


SPSla 


563.6 


546.7 


569.0 


546.6 


608.5 


97.0 


10.1.5 


1437.7 


1382.3 


1439.7 


1376.9 


1568.6 


291.3 


P 19MSSM1A, A 


339.6 


394.8 


348.3 


392.7 


414.7 


299.1 



Table 3: On-shell masses of the squarks, the gluino, and the lightest neutralino within the different SUSY 
scenarios considered. All masses are given in GeV. 



eigenstates. We use the resulting on-shell parameters directly as input for our calculation. Low scale 
soft input parameters for the pl9MSSMlA benchmark scenario are given in table ||. The physical 
spectrum is equivalcntly calculated with S0FTSUSY. For all considered benchmark scenarios we 
summarize relevant low energy physical masses in table |^. Due to non- vanishing Yukawa corrections 
implemented in S0FTSUSY the physical on-shell masses for second-generation squarks are slightly 
different from their first-generation counterparts. To simplify our numerical evaluation we set all 
second-generation squark masses to their first-generation counterparts. We checked that in the 
results this adjustment is numerically negligible. However, the general setup of our calculation is 
independent of this choice. For all considered benchmark scenarios the gluino is heavier then all 
light flavour squarks m g > rriq. Thus, these squarks decay only into electroweak gauginos and 
quarks. For SPSla and 10.1.5, right-handed squarks dominantly decay directly into the lightest 
neutralino Xi (due to its bino nature). In contrary, left-handed squarks decay dominantly into 
heavier (wino-like) neutralinos and charginos, which subsequently decay via cascades into Xii and 
only a small fraction decays directly into Xi- I n this paper we only investigate the direct decay of 
any light flavour squark into the lightest neutralino. For point pl9MSSMlA all neutralinos and 
charginos, but the lightest one, are heavier than any light-flavour squark. Thus, only the direct 
decay is allowed and all channels contribute equally to the signature under consideration. 

In table ^ we list all needed total decay widths of the squarks at LO and NLO, calculated as 
explained in section ||. NLO corrections in the total decay widths are of the order of a few percent 
for all three benchmark scenarios 3 . The total decay width of the gluino is calculated with SDECAY 
at LO and also listed in table |^. It is only used in the interference contributions of additional quark 



radiation channels discussed in section 3.2.2 



Besides physical quantities, in our calculation phase-space slicing and regulator parameters en- 



3 Numerically we found a substantial disagreement of the partial decay widths at NLO for pl9MSSMlA from the 
results obtained with SDECAY, despite the fact NLO corrections in SDECAY should also be based on the analytical 
results calculated in [46] and in appendix |^. We checked carefully, that our numerical results agree with these 
analytical results. 
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benchmarkpoint 


ul 




d L 


dR 


g 


SPSla 


r(°) 


5.361 


1.148 


5.253 


0.287 






r(°+!) 


5.357 


1.131 


5.255 


0.283 


6.849 


10.1.5 


r(°) 


12.47 


2.854 


12.46 


0.710 






r(°+!) 


12.31 


2.821 


12.30 


0.702 


10.04 


pl9MSSMlA 


r(Q) 


2.241- 1CT 3 


0.1625 


3.341 • 10~ 3 


3.917- 10- 3 






r(o+i) 


2.497- 10~ 3 


0.1621 


3.503 • 10~ 3 


3.912- 10" 3 


3.441 



Table 4: Leading order F' ' and next-to-leading order r ( - 0+1 ' total widths of light squarks and gluino for 
the considered SUSY scenarios. All widths are given in GeV. 



ter as inputs in the calculation of virtual and real NLO contributions, as explained in sections y 
and||. In the results shown in this article we set 8 S — 2-10~ 4 , Sg — 10~ 4 and m q = 10~ 15 GeV, both 
for production and decay. Numerically we checked carefully that varying their values our results 
remain unchanged on the level of individual distributions once jets are recombined (as explained 
below). We made sure that this holds for all terms of eq. (|2.7|) individually. 



5.2 Observables and kinematical cuts 

The physical signature we have in mind when calculating on-shell squark production and decay 
into the lightest neutralino is 2j + $Jt(+X). Numerical results of our calculation are presented in 
this spirit. In order to arrive at an experimentally well defined two-jet-signature we always employ 
the anti-/c<r [83] jet-clustering algorithm implemented in Fast Jet 3.0.2 [84]. Thus, we provide a 
realistic prediction on the level of partonic jets 4 . In general we use a jet radius of R = 0.4, as 
in the SUSY searches performed by the ATLAS collaboration [4]. CMS instead uses a radius of 
R = 0.5 [3]. We employ R = 0.5 in the distributions and signatures used by CMS (i.e. particularly 
the a? distribution as described below). Although we did not perform a systematic study, our 
results seem to be independent of this choice. After performing the jet clustering we sort the 
partonic jets by their transverse momentum p 1 , and in the following analysis we keep only jets with 

pJ i/2 >20GeV, 1^1 < 2.8, (5.2) 
pj. > 50 GeV, iTy^l < 3.0 (for CMS observables) . (5.3) 

Cuts of eq. ( |5.2| ) are used everywhere but in the observables used specifically by CMS (qt, as 
defined below), where cuts of eq. ( |5.3| ) are applied. 

Before showing results for the experimental signature 2j + ^?x(-|-V), we compare in section 



5.3.1 values for NLO total cross sections of squark-squark production, without decay included, with 



results obtained from Prospino 2. In section 5.3.2| we investigate the effect of NLO corrections, for 



different benchmark points, on the following differential distributions: 

• the transverse momentum of the two hardest jets pJ/ 2 > 

• the pseudorapidity of the two hardest jets 771/2, 

• the missing transverse energy J£t, 

• the effective mass m B s = ^2 pj + $t , 

i=l,2 



4 With the term partonic jets we mean that the jet-clustering-algorithm has been applied to events as produced 
from our calculation. No QCD showering or hadronization is included in the simulation. 
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• the scalar sum of the p T of all jets (visible after cuts of eq. (5.3)), Ht — Pii 

i=l,2(,3) 

• the invariant mass of the two hardest jets m- lm {jj), 

• the cosine of the angle between the two hardest jets cos0jj, which depends on the spin of the 
produced particles and therefore might help to distinguish SUSY from other BSM models [13], 

• cos = tanh {^j 2 - J , rjjj = rji — r]2, introduced in [12] as a possible observable for early spin 
determination at the LHC, 

• the olt variable, first defined in [85], where for hard real radiation events with three jets and 
pj > 50 GeV, these jets are reclustered into two pseudojets by minimizing the difference of 
the respective Ht of the two pseudojets, as explained in [86,87]. Furthermore, in all olt 
distributions we require Ht > 350 GeV as in [86] . 

Searches for sparticle production performed by ATLAS are based on p T , $t and m e g cuts; CMS 



instead uses «t to reduce SM backgrounds. In section 5.3.3 we examine NLO corrections in the 
resulting event rates after cuts. Explicitly we employ the following cuts used by ATLAS, 

pi > 130 GeV, P J 2 > 40 GeV, \ Vjl/2 1 < 2.8, ^{j 1/2 ,0 T ) > 0.4, (5.4) 
m cff > 1 TeV, gr/mrfr > 0.3, 

in their two-jet analysis [4]. Here, A</>(j'i/2, $t) denotes the angular separation between the two 
hardest jets and the direction of missing energy. Instead the CMS signal region [86] is defined as 

P J i/2 > 100 GeV, 1 < 2.5, \ Vj2 1 < 3.0, (5.5) 
H T > 350 GeV, H T j$T < 1-25, a T > 0.55, 

where JIt is calculated from Ht- 
5.3 Results 

5.3.1 Inclusive cross sections 

In table || we list inclusive LO cross sections and corresponding NLO K-factors for the three 
benchmark scenarios SPSla, 10.1.5, pl9MSSMlA, varying the LHC center of mass energy y/S = 
7, 8, 14 TeV. 

On the one hand, we list inclusive cross sections, ct^^, ~~, , and K-factors for just the production 
of squark-squark pairs summed over all flavour and chirality combinations. Next to K-factors from 
our calculation, K pp ^qqi , we also list K-factors obtained from Prospino 2, K^^^^° , where we 
adjusted Prospino 2 to use the same set of PDFs and definition of the strong coupling a s (/x/{) as 
in our calculation. Numerical differences between K pp ^qq> and K^^~}° are due to two reasons. 
First, in Prospino 2 NLO corrections to squark-squark production are calculated using an aver- 
aged squark mass. Second, the different treatment of real quark radiation (explained in section 



3.2.2 ) results in slightly smaller K-factors in our calculation. We checked numerically, excluding 
real quark radiation and putting all relevant squark masses to the same value, that inclusive NLO 
corrections from our calculation and results from Prospino 2 are in perfect agreement. 

On the other hand, we list cross section predictions for combined production and decay at LO 
<T 2°+ffi (+x) anc ^ ^ e corresponding NLO K-factors K 2 j+$ T (+x) calculated using eq. ( |2.7| ), again 



summed over all flavour and chirality combinations. Here, cuts defined in eq. (^2) are applied and 
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benchmark 


VS [TeV] 


_(o) 

pp— tqq' 


-Kpp-S-qq' 


I^Prospino 
pp—tqq' 




K 2 j + pr(+X) 




7 


1.02 pb 


1.37 


1.41 


0.37 pb 


1.41 


SPSla 


8 


1.49 pb 


1.35 


1.40 


0.53 pb 


1.40 




14 


5.31 pb 


1.28 


1.38 


1.74 pb 


1.36 




7 


0.90 fb 


1.57 


1.60 


0.45 fb 


1.61 


10.1.5 


8 


2.62 fb 


1.52 


1.54 


1.24 fb 


1.56 




14 


50.04 fb 


1.40 


1.41 


20.41 fb 


1.44 




7 


7.90 pb 


1.40 


1.37 


6.31 pb 


1.50 


P 19MSSM1A 


8 


10.48 pb 


1.39 


1.37 


8.35 pb 


1.50 




14 


29.01 pb 


1.34 


1.39 


22.60 pb 


1.47 



Table 5: LO cross sections a^ p _^-~, and NLO K-factors for inclusive squark-squark production from our 
computation, Kpp^^gi, from Prospino 2, LO cross sections of inclusive combined squark-squark 

production o^j+prC+X) an< ^ corres P° n ding K-factor K 2 j+p T (+x), f° r the three benchmark scenarios SPSla, 
10.1.5, pl9MSSMlA and center of mass energies Ss = 7,8, 14 TeV. In the combined predictions cuts of 
eq. ((5.2]) are applied. 



cross sections are strongly decreased due to the branching into the lightest neutralino. All K-factors 
of the combined process are bigger then the corresponding K-factors of inclusive production. For 
point 10.1.5 (and thus rather heavy squarks) these increasements are less then 5 %. For the other 
two benchmark points (and thus smaller squark masses) increasements in the K-factors can be of 
the order of 10 % and increase with higher center of mass energies. In general, however, K-factors 
decrease with higher center of mass energies and increase with higher masses, both, for inclusive 
production and combined production and decay. 

In table [| we compare the inclusive production with combined production and decay for bench- 
mark point 10.1.5 and a center of mass energy of \/~S = 14 TeV. Here, we list results for individual 
flavour and chirality combinations. In general, our calculation is setup to treat all 36 + c.c possible 
flavour and chirality combinations independently. However, for simplicity and to save computing 
time we always sum combinations with identical masses and matrix elements into 16 channels, both, 
for production and in the combination. This categorization follows the four possibilities discussed 



in section 3.1. For example, the channel ulUl also includes clCl (and as everywhere else in this 
paper, the charge conjugated processes). Similarly, the channel ul<1l also includes 5l§l, ul§l and 
clcIl] and the channel ulCr also includes urCl- K-factors, both, in just the production and in 
the combined result, vary by up to 15 % between different channels. Thus, an independent treat- 
ment seems adequate, as in general, squarks of different chiralities and thus different channels have 
very different decays and kinematic distributions. This can easily be seen from the very different 
order of magnitude of the different c^+^r (+X ) ^ n ^ a ^ e H ■ As already seen in table ||, K-factors 



increase comparing inclusive production and the combined result (where the cuts given in eq. (5.2) 
are applied). 

5.3.2 Differential Distributions 

Now we turn to the investigation of differential distributions in various observables. 

First, we compare the differential scale dependence between our LO and NLO prediction. We 
do this by varying at the same time renormalization and factorization scale between /j/2 < /j < 2fj,, 
with /i = Wiq. For SPSla and an energy of VS = 14 TeV resulting LO and NLO bands are shown 
in blue and red in figure |^ for differential distributions in pj,pj, 771, 772, $t and Ht- Clearly, in 



- 17 - 



channel 
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38 
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8 9^ 




1 
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40 
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±.00 


1 84 


1 


39 


UlUr 


no 




8.49 




1 49 


7 78 . 1 n~ 2 

1 . 1 O -LU 


11 ^ . 1 0~ 2 


1 


45 


dijdR 
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1 
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Uijdij 


O. ZiU 
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1 44 
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1 7R . 1 n -3 


9 fi9 • 1 O -3 

Z.UZi 1U 


1 
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1 
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1 


4.33-10" 


1 


1 ^9 
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o.uo ±u 


77Q.1 n -5 
1 . 10 lu 


1 


37 


a R^R 


4.29 • 10" 


1 


5.74-10" 


1 


1.34 


4.09 • 10" 1 


5.68 • 10 _1 


1 


39 


dh SL 


1.95 • 10" 


1 


2.75 ■ 10- 


1 


1.41 


5.16 • 10~ 5 


7.5097 • 10~ 5 


1 


46 


dRSR 


2.71 • 10- 


1 


3.87-10" 


1 


1.42 


2.59- 10- 1 


3.82 


1 


48 


u L d.R 


2.44 




3.50 




1.44 


3.16- 10" 2 


4.67- 10" 2 


1 


48 


URd L 


2.40 




3.46 




1.44 


3.87- 10" 2 


5.70- 10~ 2 


1 


48 


u L c R 


1.69-10" 


1 


2.39-10" 


1 


1.41 


2.19- 10~ 3 


3.18- 10~ 3 


1 


46 


<1lSr 


9.51 • 10- 


2 


1.39-10" 


1 


1.46 


1.52 • 10" 3 


2.29- 10~ 3 


1 


50 


sum 


50.04 


69.86 


1.40 


20.41 


29.32 


1.44 



Table 6: For the benchmark point 10.1.5 and a center of mass energy of ■s/S = 14 TeV inclusive production 



cross sections at LO a 



(0) 



., and NLO a 

11 vv 



(0+1) 



, together with the corresponding K-factors K p 



are listed 



for all different flavour and chirality channels (as explained in the text). Also listed for all channels are LO 
a 2j+$ T (+x) arm NLO °2j+|? T (+x) predictions of combined production and decay and the corresponding 
K-factor Ko 



where the cuts of eq. (5.2) are applied. All cross sections are given in fb. 



all considered distributions the scale dependence and thus the theoretical uncertainty is greatly 
reduced by our NLO calculation. At the same time, one should also note that large parts of the 
NLO bands are outside the LO bands. Still, for example in the p T distributions, in the high-pT 
tail the NLO bands move entirely inside the LO bands. 

Second, we investigate the change in the shape of distributions relevant for searches for super- 
symmetry at the LHC when going from LO to NLO. Here, we present distributions for a center of 
mass energy y/S = 14 TeV. Lower center of mass energies show qualitatively the same behaviour. 
For benchmark point SPSla plots are shown in figure [| for 10.1.5 in figure @ and for P 19MSSM1A 
in figure B. We present distributions in pj , pj, m e g, f,T (all in fb/GeV), where the ATLAS jet 



choice R — 0.4 and cuts of eq. (5.2) are applied. Also distributions in Ht (in fb/GeV) and in ctr (in 



pb) are displayed, where the CMS jet choice R = 0.5 and corresponding cuts of eq. (5.2) are applied. 
In the «t distribution, events are reclustered into two pseudojets and a cut of Ht > 350 GeV is 
applied. In the upper part of any plot we show each distribution at LO in black, NLO in red and in 
blue the LO prediction rescaled by the ratio, if NLO , between the integrated NLO and LO result. In 
the lower part of any plot we show the NLO divided by the rescaled LO • K Nho distribution. In this 
way we present corrections purely in the shape and not in the normalization of the distributions. 
For SPSla and 10.1.5 corrections are qualitatively very similar and rather flat for pj , p^ and fir- 
Corrections in the (inclusive) Ht distribution grow for larger Ht and can be sizeable. This can 
be explained from the high-p T behaviour of the contribution from hard real gluon radiation to this 
observable. Corrections to the shape of the «t observable change sign at the physical boundary [85] 
olt = 0.5 and fall off continuously in the signal region ar > 0.55. 
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Figure 5: Differential distributions in pJ,P2 (in fb/GeV), 771,772 (in pb/GeV), fix and Ht (in fb/GeV) 
for SPSla and \/S = 14 TeV, where for, both, LO (blue) and NLO (red) the common renormalization and 
factorization scale fi is varied between /_t/2 < /1 < 2/j, with /j = fng. 



Looking at the distributions of pl9MSSMlA in figure g a completely different behaviour of the 
NLO corrections cannot be missed. The tail of the pj, pj, m c g and $t distributions completely 
departs from the LO predictions. This can be understood from the following considerations. Due to 
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the small mass splitting between squarks and the x? f° r benchmark point pl9MSSMlA, jets from 
squark decays tend to be soft. Now, the p T of an additional jet (which can not be distinguished from 
the decay jets) from hard gluon radiation in the production can easily be of the same order as the 
ones from squark decays and result in the given distortions. Such a behaviour for compressed spec- 
tra was already partly discussed in [66] , where sparticle production and decay including additional 
hard jets matched to a parton shower was investigated. We verified our findings by comparing LO 
predictions plus real hard gluon radiation in the production stage with a corresponding calculation 
performed with MadGraph 5 [88]. The tail of the considered distributions can adequately only be 
described by additional gluon radiation and, corrections calculated here have effectively to be con- 
sidered as LO for these phase-space regions. Still, only our full NLO calculation allows a consistent 
treatment of the entire distributions. Turning to the olt distribution, clearly shapes of LO and 
NLO prediction are different and, here, we refrain from showing explicitly corrections in the shape 
or a rescaled LO prediction. 

Next to the distributions shown in figures |[ and ||, we also investigated pseudorapidity dis- 
tributions of the two hardest jets 771/2- Here, in the relevant region \i]i/ 2 \ < 3.0 corrections in the 
shapes are always smaller than about 5 % for all benchmark scenarios and energies. 

In figure ^ we turn our attention towards angular distributions between the two hardest jets. 
On the left we show distributions in the invariant mass of the two hardest jets mi nv (jj), on the 
right distributions in the cosine of the angle between the two hardest jets cosO^- are presented. 
Again results are shown for all three benchmark points and a center of mass energy = 14 TeV. 
Corrections in these distributions can be quite large. In general, in the full NLO results one observes 
an increase for small angles between the two hardest jets (up to 20 % in the cosQjj distributions). 
In the high-invariant-mass tail for SPSla and 10.1.5 corrections are negative and grow to 40 % in 
the considered invariant mass range. In the invariant mass distribution of pl9MSSMlA we observe 
the same deviation of the NLO result from the LO shape as already discussed above. 

Finally, in figure ^ we investigate NLO corrections to the cos O distribution for the benchmark 
points SPSla (top left), 10.1.5 (top right) and pl9MSSMlA (bottom) at a center of mass energy of 
VS = 14 TeV. Corrections up to 15 % are observable. Still, the general shape and thus the potential 
for extraction of spin information about the intermediate squarks seems to be robust under higher 
order corrections. 
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Figure 6: Differential distributions of benchmark point SPSla at a center of mass energy s/S = 14. In 
the upper part of the plots we show in black LO, in red NLO and in blue LO distributions rescaled by the 
ratio 7f NLO between the integrated NLO and LO results. In the lower part of the plots NLO corrections in 
the shapes are shown, defined as the full NLO divided by the rescaled LO • A" NLO distribution. From top 
left to bottom right we show differential distributions in pj , pj, m e g, $t, Ht (all in fb/GeV) and in cut 
(in pb). 
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Figure 7: Differential distributions of benchmark point 10.1.5 at a center of mass energy v5 = 14. In the 
upper part of the plots we show in black LO, in red NLO and in blue LO distributions rescaled by the ratio 
K NLO between the integrated NLO and LO results. In the lower part of the plots NLO corrections in the 
shapes are shown, defined as the full NLO divided by the rescaled LO • K NIj ° distribution. From top left to 
bottom right we show differential distributions in pj , pj, m e ff, $t, Ht (all in fb/GeV) and in «t (in pb). 
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Figure 8: Differential distributions of benchmark point pl9MSSMlA at a center of mass energy s/S = 14. 
In the upper part of the plots we show in black LO, in red NLO and in blue LO distributions rescaled by 
the ratio /sT NLO between the integrated NLO and LO results. In the lower part of the plots NLO corrections 
in the shapes are shown, defined as the full NLO divided by the rescaled LO • K Nh ° distribution. From top 
left to bottom right we show differential distributions in pj , pj, m e g, fir, Ht (all in fb/GeV) and in cut 
(in pb). 
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Figure 9: Distributions in the invariant mass rrn nv (jj) (in fb/GeV) and the cosine of the angle between 
the two hardest jets cosO.,., (in pb) for the benchmark points SPSla (top), 10.1.5 (middle), pl9MSSMlA 
(bottom) and a center of mass energy y/S = 14 TeV. In the upper part of the plots we show in black LO, 
in red NLO and in blue LO distributions rescaled by the ratio if NLO between the integrated NLO and LO 
results. In the lower part of the plots NLO corrections in the shapes are shown, defined as the full NLO 
divided by the rescaled LO ■ NLO distribution. 
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Figure 10: Distributions in cos O (in pb) for the benchmark points SPSla (upper left), 10.1.5 (upper right) 
pl9MSSMlA (bottom) at a center of mass energy of VS = 14 TeV. In the upper part of the plots we show 
in black LO, in red NLO and in blue LO distributions rescaled by the ratio if NLO between the integrated 
NLO and LO results. In the lower part of the plots NLO corrections in the shapes are shown, defined as 
the full NLO divided by the rescaled LO ■ if NLO distribution. 
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5.3.3 Event rates 



After investigating inclusive cross sections and differential distributions, we now proceed to event 
rates, i.e., fractionally integrated differential distributions. By this study, we want to quantify a 
possible impact of our calculation on current searches for supersymmetry and future measurements 
of event rates at the LHC. 

In table we list cross sections after applying cuts of eq. ( dA ) and in table || cross sections after 
applying cuts of eq. (5.5). We show LO and NLO cross sections for all three benchmark points and 
all three energies together with resulting K-factors. For comparison we again list inclusive K-factors 
of just production, already shown in table |^. From these results, a fully differential description of 
all squark and gluino channels including NLO effects in production and decay seems inevitable for 
a conclusive interpretation of SUSY searches (or signals) at the LHC. For any study of compressed 
spectra, like pl9MSSMlA, this seems to be eminent. Furthermore, as already suggested in [89] 



and expected from the differential distributions shown in section 5.3.2 , in particular interpretations 
based on cut seem to be highly affected by higher order corrections. 



benchmarkpoint 


Energy [TeV] 


ATLAS 


N (0+X) 

ATLAS 


-^C/Vatlas 






7 


0.066 pb 


0.083 pb 


1.26 


1.37 


SPSla 


8 


0.097 pb 


0.121 pb 


1.25 


1.35 




14 


0.347 pb 


0.424 pb 


1.22 


1.28 




7 


0.313fb 


0.503 fb 


1.61 


1.57 


10.1.5 


8 


0.861 fb 


1.344 fb 


1.56 


1.52 




14 


13.82 fb 


19.77 fb 


1.43 


1.40 




7 


0.140 fb 


20.76 fb 


- 150 


1.40 


P19MSSM1A 


8 


0.339 fb 


37.96 fb 


- 110 


1.39 




14 


0.0044 pb 


0.264 pb 


- 60 


1.34 



Table 7: LO A^las and NL0 N AThls 

cross section predictions and K-factors A"jv ATLAS for the three 
benchmark scenarios SPSla, 10.1.5, pl9MSSMlA and center of mass energies VS = 7, 8, 14 TeV where the 



cuts of eq. (5.4) are applied. For comparison we also list is the inclusive NLO production K-factor K p 



already shown in table 



benchmarkpoint 


Energy [TeV] 


JV CMS 


iv CMS 


K-Ncms 


Kp p ^.qq> 




7 


0.112 pb 


0.141 pb 


1.26 


1.37 


SPSla 


8 


0.157pb 


0.197pb 


1.25 


1.35 




14 


0.488 pb 


0.614pb 


1.26 


1.28 




7 


0.201 pb 


0.261 pb 


1.30 


1.57 


10.1.5 


8 


0.542 fb 


0.674 fb 


1.24 


1.52 




14 


8.129 fb 


8.884 fb 


1.09 


1.40 




7 


10~ 6 pb 


0.095 pb 


O(10 4 ) 


1.40 


P 19MSSM1A 


8 


10~ 6 pb 


0.151 pb 


O(10 4 ) 


1.39 




14 


2 • 10" 5 pb 


0.687 pb 


O(10 4 ) 


1.34 



Table 8: LO AT^g and NLO A^^g ; cross section predictions and K-factors Kn cms for the three benchmark 
scenarios SPSla, 10.1.5, pl9MSSMlA and center of mass energies \fS = 7,8,14 TeV where the cuts of 
eq. (5.5) are applied. For comparison we also list is the inclusive NLO production K-factor K pp ^qq' already 
shown in table || 



r(o+i) 
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6. Conclusions 

In this paper we presented a study of squark-squark production and the subsequent squark decay 
into the lightest neutralino at the LHC, providing for the first time fully differential predictions 
for the observable signature 2j + ]£ T (+X) at NLO in QCD. For the calculation, each flavour and 
chirality configuration of the squarks has been treated individually, allowing the combination of 
production cross sections and decay distributions in a consistent way. 

We studied inclusive cross sections, differential distributions for jet observables, and experi- 
mental signatures with cuts, illustrating the effect of the NLO contributions. In general, the NLO 
corrections are important, and in particular NLO effects going beyond a rescaling with a global K- 
factor cannot be neglected for setting precise limits on the sparticle masses and model parameters, 
respectively; they become specially important in model classes with compressed spectra. At the 
same time, the theoretical uncertainty on the level of differential distributions is highly reduced by 
our calculation. 

Although the present study is dedicated to the simplest squark decay mode, the fully differential 
description facilitates also the study of more complex final states including electroweak decay chains 
at the same level of the QCD perturbative order. Since the calculational framework is fully exclusive 
also with respect to flavour and chiralities, it can easily be merged with the electroweak contributions 
of LO and NLO. 
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A. Diagrams of NLO corrections 



Here, for completeness, we display all relevant diagrams used in our NLO calculation of squark- 
squark production. The contribution of some of them vanish under the assumption m q = 0. For 
example, this is the case for the 5 th diagram on the 1 st line when a ^ b; any helicity state of the 
quark in the propagator can interact either with g ija or with qj, b but not with both of them. 



Qi Qi,a Qi Qi,a Qi Qi,a qi Qi,a Qi Qi,a 




Qi Qi,a qi Qi,a qi Qi,a qi g Qi,a qi g Qi,a 




* qj Qj,b 



Figure 11: Loop diagrams contributing to all flavour and chirality structures of squark-squark production. 



- 28 - 



Qi.a 




qi 9 



- I 9 

Qi,a i 



-» * 



&,6 



Figure 12: Loop diagrams contributing only for squarks with equal flavour. 




Figure 13: Real gluon radiation diagrams contributing to all the flavour and chirality structures. Addi- 
tional diagrams contributing only for equal flavour squarks are obtained from the above ones by a simple 
crossing of the initial state quarks. 



- 29 - 



B. Formulae for soft and collinear gluon radiation 
B.l Soft radiation 



The soft-bremsstrahlung correction factor in eq. (3.10) involves the kinematical factors C~ g , 



{tt,ut,uu) 



from ( |3.9| ) and phase space integrals [43,74]. Keeping a finite quark mass m q only in the mass- 
singular terms, yields the following expressions Zy for the production process qq' qq' , where the 
labels i,j correspond to the assignment q — > 1, q' 2, q — > 3, and q' 4, 



7 — 

7ji — 
Il2 = 

Z34 = 



In 



«12 



1, (l^h 



for* = {1,2}, 
for i = {3,4}, 



(tt) 



3 99 3 qq 



(uu) 



3 ?f 3 9?' 



1 n( u ) ^ r<( uu ) L r<( ut ) 
s 



d / i=l,2 



iln 2 ^ 12 



2 1 



1 , 

6 -W 



i=3,4 
„2 



In 



l + A 



In 



/4(A£) : 



2Li 2 



2ft 
1 + 



5 la' 



l~ft 
1+ft 



In 



9 9 



A 2 



l m 2 feVlln 2 



2Li 2 1 



) - 2Li 2 



1 2 m {l 



(i-ft) 



2 \ m i / 2 
for i + j = 5 , 



1 + A 



In 



2 2 
Vfl-Vfl- 



In 



/ 4(Aff) ; 



_ 2 m 2 ( fl| ) _ A ln 2 ^ 1 - ft 



2Li 2 l -^(l + ^ -2Li 2 



for i + j = 4 or i + j = 6 , 

(B.l) 



with 



ZPi-Pj, A = Wi\/Pi, Vij = ^1 - 4m 2 m 2 /s 2 - 



For the decay process q —> q\ , the corresponding expressions for in the decay width 



eq. (4.7) read as follows, with the the label assignment q — > 1, q — > 2 (and x° ~> 3), 



In = 

Zl2 = 
I22 = 



- 2 



A 2 



7T~ 



(B.2) 



In 



4(P§) 2 



B.2 Collinear radiation 



This appendix collects the entries in eqs. ( 3.12 ) and ( 3.13 ) for hard collinear gluon radiation. They 



are taken from [75] with the replacement aQ q — > (4/3)a s . 

The (unintegrated) parton luminosity entering both the collinear gluon emission ( 3.12 ) and the 
subtraction term ( 3.13 ) is given by 



coll 



fi(~^ F ) fq'i-^Fj + fq( x ,VF) fq'(—, Vf) ■ (B.3) 

\ Z f *> iC / V XZ ' - 
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The partonic cross section for the collinear emission of a gluon into a cone around the emitting 
particle with opening angle A9 can be written as follows, 



j -i. coil-cone 
a0 qq' ^qq' 



1 + Z 2 , ( S 5, 

-j log 

1 — z 



2m 2 z 



2z 
1 - z 



(B.4) 



with Sg = 1 — cos (Ad) ~ A9 2 /2. The variable z is the ratio between the momenta of the emitter 
parton after and before the emission. 

The z-dependent part of the subtraction term at the parton level is given by 



1 + z 2 
1 - z 



log 



t' 2 F 



(1 — z) 2 m 2 



1 + z' 



1 - z 



involving the factorization scale fiF, whereas the z-independent part 



1 - log(^) - log"(S s ) + ( log(d s ) + - ) log 



(B.5) 



(B.6) 



contains also the soft-gluon phase space cut 5 S = 2AE/t/s. 

In an analogous way, the (unintegrated) parton luminosity for collinear quark radiation in 
cq. (3.18) is given by 

£ ijk (T,x,zy° n -*™ A = g g,^) [S lk fj^lip) + S jk , (B.7) 

involving also the gluon distribution function g(x,fip)- 

Collecting in one formula the collinear cone emission of a quark and the subtraction term (last 
term in the brackets of eq. ( B^ ) ) , one finds 



i-coll-quark / \ t~(0) 



s{l-z) 2 5 t 
2m 2 z 



+ 2z(l -z)- log (4) 



(B. 



with the gluon-quark splitting function P qg (z) — [z 2 + (1 — z) 2 ]/2. 
Finally, the parton luminosity for non-collinear quark radiation is 

"* dx 



Ujk{T) 4 = / — 9{X,fJ.F) 



>ik 



(B.9) 
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C. NLO corrections to the squark decay width 



In section 4.3 we need the NLO corrections to the decay width T(q — > qx), where \ is either a 
neutralino or a chargino, and q a light quark. Here we perform the derivation of the correction 
factor, following the steps of the former calculation [46], but keeping explicitly the dependence on 
the masses m q for the collinear singularities and A for the IR singularities. The NLO corrections in 



eq. ( 4.10 ) can be expressed in terms of a form factor _FQ CD , which receives four contributions 

F QCD =F + p + p ct + p r 



(C.l) 



namely loop corrections involving gluons (F g ) and gluinos (F g ), the counterterm contribution {F ct ), 
and the contribution from real gluon emission (F r ) . 

Keeping m q and A as independent mass parameters for the singular terms, F g and F ct can be 
written as follows, 



F n = 



F r , = 




1 







| + 1 - lo S 


Pi 




\ m %) 




A 1 lm\\ I A 2 



( 7 2 -l)l0g 
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"27-I 




+ 4 


i(l-7) 2+ J 



log(l - k) - log^l - k) + log(l - k) - Li 2 (k) , 




4(1-7) 



2 
2 




15 



log(7) , 



(C.2) 



(C.3) 



where k = m^/m-, 7 = m|/m|, and A denotes the UV divergence, cf. eq. (3.7). 

F g is free of soft, collinear, and UV singularities, hence it is not affected by the choice of 
regulators; it is identical to the result in [46] and we do not repeat it here. 

The part from real gluon emission, integrated over the full phase space, can be expressed with 
the help of the bremsstrahlung integrals given in [74] , evaluated in the limit m q = except for the 
mass-singular terms. The fully integrated decay width for q — > qx9 can be written as follows, 



_ r (o) 



r ' 1— /•' 

q->qX9 - 1 q^qx ' 3 n rr ' 

2 



(C.4) 



F r = 



2 2 

m| - m| 



[2{m\ - ml) (m? 7 00 + m 2 q I lx + I Q + h) - 2(m| - m?) 2 I 01 - I - I°] 



with r^ ' from eq. (4.1) and ([f.8|). The phase space integrals / = I(m,q, m q , m x ) are given by 
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Ia = 
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Am 2 - 



2 2 c\ 2 i / 
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mg + 2m^ log 
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2m| log 



16 TO? 
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(C.5) 



With these expressions, eq. flC.4|) yields J 71 ,, for the real gluon part of _fQ cd in eq. (p.l 



-5 + 8k - 3k 2 - 8k bg(«) + 6k 2 log(/c) 



8(1 -k) 2 



4-^-|log(l 



1 



log(l-K) 2 -log 2 (^)+21og(^) 



- log ^ j + 2 log(l - k) log (^i) + log 2 (^i) + Li 2 (1 - k) 



(C.6) 



Combining all four contributions in eq. (C.l) we obtain a compact analytical expression for the 
form factor f^ CD , which agrees with the result in [46], 



f qcd 



1 /4 7 2 - 277 + 25 3k-5 
8 I 



7-1 

3 7 2 -4 7 + 2 



4(1-7) 



2 log (7) 
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2U 9 (k) 
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(C.l) 



where for K7 < 1 the function / is given by 



/ = Li, 



7-1 
7K — 1 



Li 2 ( k— — 7] - Li 2 ( 7 + K - 2 ) + Li 2 f/c' 1 ' 



7K — 1 



7K — 1 



7K — 1 



and for K7 > 1 one has 



-Li, 



7K - 



1 



Li 2 



7« 



7 + K-2 



Li 2 



7K — 1 
«(7~1) 



log(«) log 



7 + K-2 



— Li 2 



7K — 1 



k( 7 + k-2) 
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